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I . INTRODUCTION 


The  energetics  of  the  propellant  are  of  fundamental  importance  in 
ballistics  technology.  The  amount  of  energy  which  can  be  utilized  to 
impart  thrust  to  the  projectile  depends  upon  the  reaction  mechanism  of 
the  particular  propellant.  As  the  propellant  is  activated,  its  defla- 
gration proceeds  along  certain  reaction  paths  changing  from  one  energetic 
state  to  another.  The  energy  obtainable  from  a state-to-state  transition, 
as  well  as  the  time  a given  transition  takes,  may  be  calculated  by 
quantum  chemical  means.  In  order  to  determine  the  energy  one  may  use  the 
"transition  state"  procedure  which  relates  total  energy  differences  to 
differences  in  electronic  orbital  energies  (single-electron  eigenvalues). 
This  procedure  is  valid  for  any  states  which  are  describable  in  a single- 
electron picture;  that  is,  when  the  effects  upon  a single  electron  are, 
to  a good  approximation,  due  to  a potential  field  which  is  an  average 
field  generated  by  all  the  other  electrons  taken  in  their  entirety. 


In  this  paper  we  describe  a method  for  calculating  single-electron 

1 2 

eigenvalues  by  the  X-a  method  with  non-muffin-tin  (NMT)  corrections  ’ . 

We  calculate  and  compare  ionization  energies  using  the  X-a  transition 
state  method  (referred  to  as  NOT  or  NOT-transition)  and  the  method  of 
total  energy  differences  (referred  to  as  ASCF) , both  with  NOT  corrections, 
for  the  molecules  CO,  NO,  SN,  and  H2CO.  The  results  for  the  NO  and  SN 

molecules  using  the  method  of  total  energy  differences  have  been  previous- 
ly reported3.  Substantial  improvement  is  obtained  over  muffin-tin  (OT) 
calculations  for  the  NOT  X-a  total  energy  calculations.  The  results 
compare  very  well  with  experiment  and  with  other  methods  which  have 
been  advanced  for  calculating  X-a  transition-state  energies  such  as 

4 

the  overlapping-spheres  (OS)  method  , the  discrete-variational  method 

(DVM)3,  and  the  method  of  Sambe  and  Felton6.  The  present  NOT  X-a  method 
is  comparatively  economical  and  can  be  extended  to  larger  molecules,  just 
as  was  the  original  multiple-scattering  X-a(MSX-a)  method  in  the  muffin- 
tin  approximation  (referred  to  OT-transition) . The  work  reported  here 
was  done  primarily  to  test  the  method  and  to  determine  its  usefulness. 


J.B.  Danese  and  J.W.D.  Connolly,  J.  Chem.  Phys.  61,  3063  (1974); 
J.B.  Danese,  J.  Chem.  Phys.  61,  3071  (1974). 

J.B.  Danese,  Chem.  Phys.  Lett.,  45,  150  (1974). 

D.A.  Ringers,  "Electronic  Excitation  in  NO  and  SN  by  the  Non- 
Muffin-Tin  X-Alpha  Method,"  USA  Ballistic  Research  Laboratory, 
Report  No.  BRL  R 1968,  March  1977.  (AD  #A037802) 

N.  Rosch,  W.G.  Klemperer,  and  K.H.  Johnson,  Chem.  Phys.  Lett. 

23,  149  (1973). 

D.E.  Ellis  and  G.S.  Painter,  Phys.  Rev.  B2,  2887  (1970);  E.J. 
Barrends,  D.E.  Ellis,  and  P.  Ros,  Chem.  Phys.  2,  52  (1973). 

H.  Sambe  and  R.H.  Felton,  J.  Chem.  Phys.  62,  1122  (1975). 
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I I . THEORY 


The  relationship  between  the  MT  approximation  and  the  NMT  corrections 

for  the  X-a  total  energy  and  eigenvalues  has  been  described  elsewhere*. 

In  essence,  the  correction  to  the  eigenvalues  is  obtained  as  a first- 
order  perturbation  correction  to  the  MT  one-electron  eigenvalues.  The 
perturbation  is  the  difference  between  the  total  one-electron  potential 
obtained  using  the  orbitals  of  the  MT  one-electron  equation  [the  set 
(Q) ] and  the  MT  average  of  the  potential  obtained  using  the  MT  density. 

In  what  follows  a bar  (-)  over  a quantity  will  indicate  a MT  function 
and  a A preceding  a symbol  will  indicate  a NMT  function.  The  X-a  one- 
electron  differential  equation  for  the  orbital  u^  of  the  eigenvalue 

is 


[-V2  + v ( [p ] , r)]  ut(r)  = eiui(r)  . (1) 

The  usual  MT  X-a  one  electron  equation  is 

[ -V2  + v ([  p]  , r)]  TkC?)  = e^.  (r).  (2) 

Note  the  different  functional  dependences  upon  p in  the  two  equations  (p 
versus  p).  As  is  shown  in  reference  1,  the  MT  one-electron  potential  is  a MT 
function  and  also  depends  functionally  on  p.  The  X-a  potential  consists 
of  the  three  types  of  contributions 

vflpj.r)  = vNe(r)  + vc  ([  p]  , r)  + j caP  (r) 1/3,  (3) 

where  the  first  term  is  the  nuclear-electron  potential,  the  second  is 
the  electron-electron  Coulomb  interaction  potential,  and  the  last  is 
the  X-a  approximation  to  the  exchange  potential.  Taking  the  difference 
of  the  two  potentials  of  Eq.  (1)  and  Eq.  (2),  we  get  for  the  pertur- 
bation 


Av(r)  = AvNe  (r)  + Avjr)  + Avexc(r).  (4) 

where  the  first  two  terms  are  the  NMT  analogues  of  the  first  two  terms  in 
Eq.  (3),  and  Av£xc  is  given 

Avexc(?)  = Ica(p(^)1/3  -~r)1/3),  (5) 

(here  c^*  -9a(3/8ir)  *^3  and  a is  the  exchange  parameter2).  The  definition 
of  the  Coulomb  potential  of  Eq.  (1)  is,  in  rydbergs. 


7. 


K.  Schwarz,  Phys.  Rev.  BS,  2466  (1972). 


Mtl  dp 
I?-*  I 


(6) 


vc([  p]  ,r)  = 


We  may  break  up  p and  into  their  MT  and  NMT  parts, 
p(r)  =T(r)  + Ap  (r) 

and 

vc([p],r)  = vc(M  ,r)  + vc([~],r)  + ”c([Ap],r)  + vc([Ap],r).  (7) 

Hence  Av£  of  Eq.  (4)  is 

Avc(r)  = Avc(fp],r)  + T([Ap],r)  + Avc([Ap],r).  (8) 

With  this  definition  the  first  order  perturbation  Ae.  to  the  MT  eigenvalue 
is  given  by  1 

Aei  ^/*{Ap i(r)AvNe(r)  + Api(r)  Avc(f7]  r)  +TiCr)v'c([  Ap]  ,r) 


+ Ap.(r)Avc([  Ap] ,r)  + pA(r)  Avexc 


(r)  | 


dr. 


(9) 


Here  use  has  been  made  of  the  fact  that  the  integral  of  a MT  function 
with  a NMT  function  is  always  zero(l).  (The  exchange  part  is  not  of 
this  form  as  may  be  seen  by  examining  the  definition  of  Avexc>  Eq.  (5). 

It  is,  however,  a three-dimensional  integral.)  The  third  term  of  Eq. 

(9)  may  be  rearranged  in  a different  and  computationally  more  convenient 
way. 


([  Ap]  ,r)dr 


-/ 


Ap(r)Avc([pi]  r)  dr. 


(10) 


With  this  rearrangement  all  the  integrals  of  Eq.  (9)  are  three  dimen- 
sional except  the  last  one.  Being  three  dimensional,  they  are  amena- 
ble to  evaluation  using  three  dimensional  Gaussian  integration  as 
described  in  reference  2.  The  last  integral  is  basically  a six  dimen- 
sional integral,  which  for  a general  molecule  cannot  be  reduced  further, 
(although  certain  simplifications  do  occur  for  highly  symmetric  molecules). 
As  part  of  this  study  we  will  determine  the  importance  of  this  term  in 


The  ionization  energies  are  calculated  here  using  the  well-known 

8 9 

transition  state  formulated  for  the  X-alpha  theory  by  Slater,  et  al.  ’ , 
in  which  the  occupation  of  the  level  in  question  is  reduced  by  one-half 


8.  J.C.  Slater  and  J.H.  Wood,  Intern.  J.  Quan.  Chem.  48,  3 (1971); 
J.C.  Slater  and  K.H.  Johnson,  Phys.  Rev.  B5,  844  (1972). 

9.  J.C.  Slater,  Quantum  Theory  of  Molecules  and  Solids,  Vol.  4 
(New  York:  McGraw-Hill  Book  Co.,  Inc.,  1973). 
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electron.  This  approach  allows  the  calculation  of  total  energy  differ- 
ences by  an  approximation  involving  only  single-electron  orbital  energies 
(eigenvalues);  i.e., 

E (n.)  - E (n.  .)  =»  e.(n.  - 1/2),  (11) 

01  e l-l  ii 

where  E is  the  total  electronic  energy  of  a system  having  an  initial 

electron  in  the  il  orbital.  An  electron  is  excited  out  of  this  i 
orbital  into  the  continuum,  and  Eg  is  the  total  electronic  energy  of 

this  final  state.  The  single-electron  eigenvalue,  e^,  which  is  cal- 
culated with  the  occupation  of  the  i1^  orbital  halfway  between  that  of 
the  original  and  final  states,  is  a good  approximation  to  this  total 
energy  difference,  especially  forjgxcitations  involving  valence  electrons. 
(For  core  electrons  another  term  appearing  on  the  right-hand  side 
of  Eq.  (11)  may  be  significant.)  The  advantage  of  this  procedure  is 
that  single-electron  orbital  energies  are  generally  more  readily  cal- 
culable, and  therefore  less  costly,  than  total  energy  differences. 

For  an  optical  transition  one-half  electron  is  transferred  from  the 
lower  to  the  upper  level.  There  is  no  problem  in  combining  this 
concept  with  the  NMT  correction  described  above:  the  effect  of  the 
change  in  occupation  numbers  enters  through  the  change  in  the  density 
and  therefore  in  the  potential  terms  in  Eq.  (9).  For  an  optical 
transition  it  is,  of  course,  necessary  to  calculate  the  NMT  correction 
to  both  levels  involved. 


III.  RESULTS 

The  multiple  scattering  X-a  parameters,  internuclear  separations, 
and  alpha  exchange  parameters  used  in  the  calculations  reported  here 
are  given  in  Table  I.  Touching  spheres,  whose  radii  were  chosen  to  be 
proportional  to  their  Slater  atomic  radii,  were  used  for  the  atomic 
regions.  The  a values  were  taken  from  atomic  calculations  for  the 
atomic  regions,  with  an  averaged  value  for  the  intersphere  and  outer 
sphere  regions. 


10.  D.A.  Ringers,  "Xa  Calculation  of  Transition  Energies  in  Multiply 
Ionized  Atoms,"  USA  Ballistic  Research  Laboratory  Report 
BRL  MR  2766,  July  1977.  (AD  #A043091) 
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TABLE  I. 


Calculational  Parameters. 


Molecule 

r 

e 

°N 

“o 

“c  “s 

“IS 

NO 

2.1758 

.74725 

.74118 

.74422 

CO 

2.17 

.74118 

.75331 

.74725 

SN 

2.8765 

.74952 

.72183 

.73274 

h2co 

2.2825  rc 

-0 

.74118 

.75331 

.75006 

2.1090  r 

-H 

< OCH  133.3° 

In  Tables  II  through  V and  Figures  1 through  4 we  give  the 
ionization  energies  as  determined  by  A SCF  and  by  transition  state 


TABLE  II.  CO  Ionization  Energies  (eV) . 


Level 

MT-Transition 

A SCF11 

NMT-Transition 

Experiment12 

0 Is 

551.1 

550.8 

539.8 

542.1 

c Is 

304.9 

303.4 

302.9 

295.9 

3 a 

32.7 

36.3 

35.5 

38.3 

4 a 

19.1 

20.5 

19.4 

20.1 

1 TT 

18.1 

17.1 

16.3 

17.2 

5 a 

12.9 

14.7 

13.9 

14.5 

11.  D. 

.A.  Ringers  and  J.B. 

Danese, 

"Quantum  Chemistry, 

Chemisorption, 

and  Erosion.  I.  Properties  of  Diatomic  Molecules",  USA  Ballistic 
Research  Laboratory  Report  BRL  R 2023,  October  1977.  (AD  #A050249) 
12.  K.  Siegbahn,  C.  Nordling,  G.  Johansson,  J.  Hedman,  P.F.  Heden, 

K.  Hamrin,  U.  Gelius,  T.  Bergmark,  L.O.  Werme,  R.  Manne,  and 
Y.  Baer,  ESCA  Applied  to  Free  Molecules  (Amsterdam:  North-Holland 
Publishing  Co.,  1971). 


r 


Figure  2.  Ionization  energies  for  NO. 
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TABLE  III.  NO  Ionization  Energies  (eV) . 


1 


Experiment12 

A SCF11 

Transition 

state  (non-spin-polarized) 

Level 

Level 

NMT 

mt 

0 1st 

544.0 

540.7 

0 Is 

543.7 

551.6 

0 1st 

543.3 

540.5 

N 1st 

411.8 

409.6 

N Is 

415.2 

419.3 

N 1st 

410.3 

408.3 

3 at 

43.8 

38.5 

3 a 

37.4 

34.4 

3 at 

40.6 

37.7 

4 at 

23.3 

23.1 

4 a 

21.7 

22.0 

4 at 

21.7 

22.0 

5 at 

18.5 

18.8 

5 a 

17.7 

14.9 

5 at 

16.7 

17.7 

1 Tit 

14(aV)a 

18.4 

1 TT 

17.2 

18.8 

1 7ft 

17.4 

2 nt 

10 

10.1 

2 7T 

8.4 

10.3 

This  is  the  highest  lying  state  of 

1 7t  orbital.  Values  for  the  *A,  , 

given. 

the 

V 

five  multiplets  arising  from  the 

3 

, and  A multiplets  are  not 

TABLE 

IV.  SN  Ionization  Energies  (eV) . 

Level 

A SCF11 

Transition 

Level 

state 

(non-spin- 

NMT 

-polarized) 

MT 

5 at 

28.7 

5 a 

28.6 

26.2 

5 at 

27.9 

6 at 

19.0 

6 a 

18.0 

18.8 

6 at 

18.2 

7 at 

14.6 

7 a 

13.7 

12.6 

7 at 

14.1 

2 ift 

13.6 

2 7T 

12.8 

14.3 

2 nt 

12.7 

3 nt 

8.6 

3 TT 

7.1 

8.9 

14 


TABLE  V. 

H^CO  Ionization  Energies 

(eV). 

Level 

Experiment^  3 

A SCF 

Transition 

NMT 

state 

MT 

0 Is 

542.1 

546.3 

529.3 

547.9 

C Is 

300.3 

296.7 

305.4 

1 al 

33.1 

32.8 

31.4 

1 bl 

10.8 

16.8 

16.5 

17.5 

1 b2 

14.4 

14.3 

13.0 

16.5 

2al 

19.9 

18.5 

20.8 

3al 

16.0 

15.8 

14.1 

17.3 

2bl 

10.9 

10.3 

9.7 

12.8 

calculations  for  the  molecules  CO, 

NO,  SN, 

and  H2C0.  In 

each  case  all 

valence  levels  are  reported  for  those  orbitals  which  are  occupied  in 
the  ground  state.  Except  for  SN,  the  core  levels  are  reported  as  well. 
The  levels  denoted  by  arrows  ( + or  + ) were  calculated  in  a spin- 
polarized  manner  using  different  orbitals  for  different  spins.  The 
Tables  also  list  the  MT-transition-state  results  so  as  to  show  the 
improvement  obtainable  with  the  NMT  procedure.  Tables  VI  and  VII  list 
the  detailed  contributions  to  the  corrections  for  the  CO  molecule. 


TABLE  VI. 

NMT  Contributions 

to  Ionization  Energies 

in  CO  (eV) . 

Level 

-ael 

aeq 

-AEt  (NMT) 

0 Is 

35.0 

5.3 

29.7 

C Is 

34.6 

3.  7 

30.9 

3 a 

29.3 

3.4 

25.9 

4 a 

32.4 

4.1 

28.3 

1 TT 

34.1 

3.7 

30.4 

5 a 

32.6 

4.6 

28.0 

13. 


D. "Turner,  C.  Baker,  A. D.  Baker,  and  C.  Brandle,  Molecular 
Photoelectron  Spectroscopy,  (New  York:  Wiley  Interscience, 


1970) . 
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TABLE  VII.  NMT  Contributions  to  Transition  State  Ionization  Energies 
in  CO  (eV) . 


Level 

-e.(MT) 

AEi(L) 

-ei(NMT) 

-e^ (Experiment) ^ 

0 Is 

551.1 

11.0 

540.0 

542.1 

C Is 

304.9 

2.1 

302.8 

295.9 

3 o 

32.7 

- 2.8 

35.5 

38.3 

4 a 

19.1 

- 0.3 

19.4 

20.1 

1 7T 

18.1 

1.8 

16.3 

17.2 

5 a 

12.9 

- 1.0 

13.9 

14.5 

Examination  of  the  results  shows  the  improvement  achieved  with  the 
NMT  linear  correction  and  the  agreement  with  experiment  in  the  ioniza- 
tion energies. 

One  of  the  purposes  of  this  study  was  to  determine  the  importance 
of  the  quadratic  density  correction  of  Eq.  (9) , which  may  be  related 
by  the  transition-state  concept  to  a difference  of  quadratic  corrections 
to  the  total  energy.  Examining  Table  VI,  we  see  that  the  quadratic 
term  for  the  total  energy  correction  has  an  average  value  of  4.1  eV 
with  the  largest  deviation  for  a valence  state  being  0.7  eV  from  this 
average  value.  Hence,  if  we  use  an  average  value,  we  make  at  most  an 
error  of  0.7  eV.  Table  VII  also  gives  the  directly  calculated  linear 
correction  to  the  eigenvalues  taken  from  Eq.  (9).  For  the  total  energy 
differences  it  will,  in  many  cases,  not  be  necessary  to  calculate  the 
quadratic  term  explicitly  for  each  total  energy  state,  or  even  for  any 
total  energy  state,  since  the  near  equality  of  this  term  for  most  total 
energy  states  involved  in  ionization  will  cause  large  amounts  of 
cancellation.  For  the  calculation  of  excitations  by  the  transition- 
state  concept,  the  quadratic  correction  to  the  single-electron  energy 
levels  will  in  general  not  be  necessary,  since  the  near  equality  of 
the  results  with  only  the  linear  correction  to  the  results  obtained 
from  a total  difference  calculation  implies  that  there  must  be  con- 
siderable cancellation  among  the  quadratic  corrections  for  single- 
electron eigenvalues  as  well  as  for  total  energy  differences.  The 
linear  correction  alone  gives  substantial  improvement  over  the  MT 
eigenvalues,  as  may  be  seen  in  the  Tables  and  Figure  1.  We  conclude, 
then,  that  for  the  identification  of  ionization  levels  the  quadratic 
term  is  necessary  only  for  levels  which  are  very  close  together.  An 
example  of  just  how  close  together  the  levels  may  be  can  be  seen  in 
the  cases  of  SN  and  NO.  The  result  using  only  the  linear  correction 
term  gives  the  correct  ordering  for  the  levels  which  are  found  to  be 


It 


separated  by  only  .025  rydberg  in  the  more  complete  calculation  (one 
including  quadratic  terms)  and  in  fact  differs  from  the  more  complete 
calculation  by  only  .006  rydberg. 

The  linear-correction  result  compares  favorably  with  the  discrete 

variational  method  (DVM)5,  Felton-Sambe6 , Hartree-Foch14  using  Koopman's 

theorem,  and  overlapping  sphere^  results  for  NO.  For  SN,  however,  the 
overlapping  sphere  results  give  an  inverted  ordering  which  has  the  7a 
level  appearing  above  the  2ir  level.  The  ordering  of  the  SN  and  NO  ioni- 
zation levels  was  found  to  be  independent  of  whether  the  calculations 
were  done  spin-polarized  or  non- spin-polarized  and  whether  the  energy 
differences  were  calculated  as  differences  in  total  energies  or  were 
approximated  by  the  transition  state  technique. 

The  core  levels  of  CO  compare  favorably  with  experiment  but  not 
so  well  as  the  valence  levels.  Further  study  is  required  to  determine 
whether  the  core  levels  might  not  be  improved  or  whether  systematics 
exist  in  the  behavior  of  these  levels  in  similar  systems. 


IV.  CONCLUSIONS 

The  inclusion  of  NMT  corrections  to  the  Xa  procedure  generally 
used  to  calculate  electronic  energy  levels  has  significantly  increased 
the  accuracy  of  the  calculated  values.  When  combined  with  the  tran- 
sition state  method,  the  procedure  is  a viable  one  for  the  calcula- 
tion of  approximate  total  energy  differences  such  as  excitation 
energies . 

It  has  been  previously  shown  that  the  quadratic  correction  term 
in  the  expression  for  the  total  energy,  while  the  most  time  consuming 
to  calculate,  is  in  general  small  and  may  usually  be  neglected  when 
calculating  total  energy  differences.  The  excitation  and  ionization 
energies  calculated  using  the  linear  correction  to  single  electron 
eigenvalues  (by  the  transition  state  method)  agree  well  with  the  values 
obtained  from  total  energy  differences.  This  implies  that  the  quadratic 
correction  for  single-electron  eigenvalues  must  likewise  be  small  and 
may  often  be  regarded  as  negligible. 

The  calculated  excitation  energies  are  somewhat  less  accurate  than 
those  determined  from  Xa  non-muf fin-tin  total  energy  differences.  The 
determination  of  an  excitation  energy  by  total  energy  differences, 
however,  requires  two  quite  accurate  calculations;  one  for  the  initial 
state  and  one  for  the  final  state.  Thus  the  slight  decrease  in 
accuracy  in  the  present  method  is  offset  by  the  need  for  only  one 
calculation  when  using  this  procedure  for  calculating  excitation 
energies . 


14.  W.M.  Huo,  J.  Chem.  Phys.  4^5,  624  (1965). 
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